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OXIDATION OF CARBON AT HIGH TEMPERATURES: 

REACTION-RATE CONTROL OR TRANSPORT CONTROL 

By Howard G. Maahs 
Langley Research Center 

SUMMARY 

In spite of a general acceptance in making predictions of ablation rate that the oxi- 
dation rate of carbon is limited at temperatures of about 1650° K or above by the rate of 
oxygen transport to the carbon surface, a large body of existing experimental data on the 
carbon-oxygen reaction indicates that this may not be so. Furthermore, these data indi- 
cate a distinct departure from simple Arrhenius behavior and even a maximum in the rate 
curve as a function of temperature. These data have been analyzed and compared with 
theoretical maximum stagnation -point transport-limited rates. This comparison has 
shown that even at very high temperatures, in many situations the rate of chemical reac- 
tion is considerably below the theoretical maximum transport -limited rate. Hence, the 
carbon-oxygen reaction may be chemically controlled in many situations at much higher 
temperatures than previously suspected. 

INTRODUCTION 

There appears to be a general, unquestioned acceptance among investigators con- 
cerned with carbon or graphite as a material for high-temperature-reentry applications 
that the oxidation rate of carbon is limited at temperatures of about 1650° K or above by 
the rate of oxygen transport to the carbon surface. (See ref. 1.) While the concept of 
transport control of a gas-solid reaction is certainly theoretically sound, the assumption 
that all carbon -oxygen reactions (or even most of those of practical interest) are inevi- 
tably transport controlled above 1650° K is inconsistent with a considerable body of 
experimental data indicating the contrary. (See refs. 2 to 22.) For a situation of trans- 
port control to exist, it is a fundamental requirement that the chemical -reaction rate far 
exceed the oxygen -transport rate to the solid-carbon surface. Analysis of the previously 
cited experimental data suggests that often this is not true for many high-temperature 
situations of practical importance. 

The purpose of this paper is to present a new analysis of the high-temperature oxi- 
dation rate of carbon for use as a material for reentry applications. The foundations of 
the present, generally accepted theory of the transport-controlled oxidation rate of carbon 



at high temperatures will be briefly discussed, followed by an analysis of the previously 
cited experimental data. Finally, it will be shown that in many situations even at very 
high temperatures, the rate of chemical reaction does not even approach the theoretical 
maximum transport -limited rate, thereby demonstrating that the carbon -oxygen reaction 
may be chemically controlled in many situations at much higher temperatures than pre- 
viously suspected. 


SYMBOLS 


G constant 

m mass-loss or reaction rate, g/cm 2 sec 

mp theoretical maximum stagnation -point transport -limited mass-loss rate 

(eq. (7)), g/cm 2 sec 

p air pressure, atm (1 atm = 101.3 kN/m 2 ) 

cl 

p^ oxygen pressure, atm 

°2 

p^. g stagnation pressure of air adjacent to the surface, atm 

R gas constant, 1.987 cal/mole °K 

r n nose radius, cm 

s linear recession rate due to reaction, cm/sec 

T absolute temperature, °K 

CURRENTLY ACCEPTED TRANSPORT-CONTROL THEORY 

The present, generally accepted theory of transport control of the carbon-oxygen 
reaction at high temperatures is largely due to Scala (ref. 23). Two key features of the 
theory are the rate at which oxygen can be transported to the surface and the rate at 
which the chemical reaction is proceeding at the surface. To obtain an estimate of the 
chemical reaction rate at high temperatures, Scala reviewed the kinetic data and reaction- 
rate expressions presented by several authors. (See refs. 24 to 28.) Pertinent informa- 
tion from Scala's data sources is presented in table 1; reaction temperature and pressure 
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TABLE 1.- SUPPORTING DATA FOR KINETIC ’’BOUNDS" OF SCALA 


Author Reference 


Range of supporting kinetic data Reactant 
T, o K p 02 , atm 


gas 


Carbon studied 


Kinetic expression presented by 


Author 


Scala (ref. 23) 


Footnote 


Gulbransen 24, 29 698 to 848 0.002 to 0.13 O 2 Spectroscopic grade; 

and Andrew Ash, 0.02%; 

Density, 1.6 to 1.7 g/cm 3 

Blyholder 25 873 to 1073 0.0000013 to 0.0002 O 2 Spectroscopic grade; 

and Eyring 


Trace impurities; 
Density, 1.3 g/cm 3 


m = 1.19 X io 2 e' 3670 °/ RT m = 1.19 X 10 2 e“ 36 700/RT ( a ) 

m = CPo 2 1/2 m = 1.23 x 10 2 e‘ 42 300 / RT (b) 


Parker 
and Hottel 


26 973 to 1373 0.002 to 0.20 Air Brush carbon; 


Frank- 27 (d) 

Kamenetskii 


(d) 


(d) 


Ash, 0.17%; 

Density, 1.55 g/cm 3 


(d) 


. 9.55 X 10 6 „-44 000/RT„ 
P 0 2 


m = — Pn„ A = 9-55 X 1 0 6 e -44000/RT (c) 


rpl/2 


41 100 , 


T-(T-Tq) MO^(T-To) 

A = Ce RT o 2 p 02 m ril = M^! e ^7 Po 2 1/2 (d) 


where where 

T 0 = 1000° K; and 1/2 s m < 2/3 


T 0 = 1000° K 


17 720 / 


Vulis 


28 733 to 993 


^ 0.21 


Air Welding carbon; 

Ash, 3.2%; 

Density, 1.5 g/cm 3 


s = 5 X loV 17 50 °/ RT m = 1^3 e R V* Vp Q (e) 


L.L) 

T* Tj. 


where 


17 500/ 1 1 


s = 31. 6e 


R IT* Tj 


T* = 1240° K 


where 


T* = 1240° K 


CO 


a The activation energy was obtained from reference 24 and the pre exponential from reference 29; both are valid only at an oxygen pressure of 0.10 atm. In reference 24 the 
order of the reaction was shown to vary with pressure, being nearly first order above 0.13 atm, zero order between 0.002 atm and 0.02 atm, and intermediate order between 0.02 atm 
and 0.13 atm. 

b Blyholder and Eyring state that their data can be expressed with a constant activation energy of 42 300 cal/mole and a reaction order of one-half, but do not present a pre- 
exponential term. They further state that above 1073° K both the activation energy and the reaction order change. Scala' s expression can be obtained by assuming an activation 
energy of 42 300 cal/mole and a zero-order reaction and by calculating a pre exponential from data given in Blyholder and Eyring's paper at an oxygen pressure of 0.000034 atm. 

c Parker and Hottel' s data were obtained at low flow velocities where the reaction rate was dependent on flow velocity and, hence, was transport influenced. Accordingly, they 
carefully specified p 0 ^ in their kinetic expression to be the true oxygen pressure at the carbon surface. The carbon studied is described in reference 30. 

d Frank-Kamenetskii’s rate expression is based on a reanalysis of Parker and Hottel’s data (ref. 26). Frank -Kamenetskii concluded that the order of the reaction lay between 
one-half and two-thirds, with an activation energy of 41 100 cal/mole; however, no preexponential term was presented. Scala's expression can be reconstructed by assuming a reac- 
tion order of one-half and an activation energy of 44 000 cal/mole and by using a parameter from Frank -Kamenetskii's paper along with figure 6 of Parker and Hottel's paper to obtain 
a preexponential term. Use of this parameter accounts for the T*/ 2 appearing in the denominator of Scala's expression. It is interesting to note that the unconventional form of 
Frank-Kamenetskii's rate expression obscures the true meaning of the activation energy in the conventional sense. 

e Although Vulis conducted his experiments in air (test conditions are given in ref. 31), the test pressure is not stated and can only be assumed to be atmospheric. Vulis pre- 
sented two rate expressions (which are virtually equivalent) but failed to specify a reaction order. The present author has not been able to reconstruct the rate expression presented 
by Scala, and the reason for its differing so markedly from those presented by Vulis is unknown. The high recession rates which are predicted by Vulis' expression, as well as the 
low reported activation energy of 17 500 cal/mole, can be attributed to catalytic accelerating effects from the very high concentration of impurities (3.2 percent ash) in Vulis' carbon 
samples. (See ref. 32.) 



ranges, reactant gas, and carbon type, as well as the associated kinetic expressions 
are also included. It is observed that the kinetic expressions attributed by Scala to his 
sources differ in some instances from those which can actually be found in those sources. 
In order to permit engineering estimates of reaction rates from this diversity of kinetic 
expressions, it was necessary for Scala to arbitrarily select a kinetic form, a reaction 
order, a preexponential factor, and an activation energy which best bracketed the data. 
The resulting kinetic bounds (in the notation of the present report) were given by Scala 
in reference 23 as 


m slow = 2. 1 8xl 0 4a-^300/RT po2l /2 

(1) 

m fast = 3.29 X 10 e e' 4400 °/ RT p 4 / 2 

(2) 


which span a reaction-rate range of four orders of magnitude. Two major facts are to be 
noted regarding these kinetic expressions: first, they are admittedly arbitrary (ref. 23); 
and second, they are all based on data obtained below 1373° K. 

Faced with the problem of estimating carbon -oxidation rates at higher tempera- 
tures, it has been customary for the engineer to extrapolate equations (1) and (2) linearly 
without regard to the temperature and pressure -range limitations of the data or to several 
cautionary statements concerning changing reaction order and activation energy given by 
the authors of the data. (See table 1 and refs. 24 to 32.) Although such extrapolation is 
convenient and often even necessary in the absence of more appropriate data, such extrap- 
olation is without theoretical justification, particularly in the case of heterogeneous reac- 
tions. (See, for instance, refs. 11, 33, and others.) Rosner and Allendorf (refs. 16 and 18) 
have shown that linear extrapolation of simple Arrhenius expressions such as equations (1) 
and (2) will, at high temperatures, often lead to predicted oxidation-rate probabilities 
exceeding those theoretically possible (i.e., predictions of more than one carbon atom 
being removed from the surface per one oxygen atom colliding with the surface). 

The reason a linear Arrhenius expression cannot be expected to adequately describe 
a gas-solid reaction is that there may be a large number of separate, competing kinetic 
steps involved (such as adsorption of the reactant, surface migration of the reactant to an 
active site, chemical combination with a carbon atom at that site, desorption of the prod- 
uct, etc.). Each of these steps may be of greater or lesser importance as the temperature 
or pressure is changed. Additional complicating kinetic steps are introduced if the reac- 
tant dissociates (as might oxygen, O 2 ), if multiple products are formed (as, for instance, 
carbon dioxide, CO 2 , and carbon monoxide, Co), if an inert diluent in the reactant gas 
competes for active sites (as might nitrogen, N 2 , in air), if there is more than one type 
of active site, and so forth. Since these kinetic steps are, in general, temperature and 
pressure dependent, both the apparent activation energy and reaction order of the 
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reaction may change as the temperature and pressure change. Therefore, although it is 
often possible, over a sufficiently narrow range of temperature and pressure, to corre- 
late heterogeneous kinetic data by a simple Arrhenius expression, such an expression 
cannot be applied with reliability outside that range unless it is somehow known that the 
predominating kinetic mechanism is the same. 

KINETIC EXPRESSIONS AT HIGH TEMPERATURES 


A number of kinetic expressions have been developed (see refs. 2 and 8 to 11) from 
high-temperature kinetic data which make an extrapolation of equations (1) and (2) unnec- 
essary. These kinetic expressions are based on data obtained at temperatures as high as 
2773° K, but at experimental conditions such that reaction-rate limitations due to oxygen 
transport were absent. One such kinetic expression is that of Blyholder, Binford, and 
Eyring (ref. 2). It is based on low-pressure kinetic data taken by the hot-carbon-filament 
technique drawn from references 3 to 7. Figure 1 shows these data as adjusted by 
Blyholder et al. (ref. 2) to a common pressure of 10”^ torr oxygen (1 torr = 133.3 N/m^). 


TEMPERATURE, °K 

2500 2000 1667 1429 1250 1111 1000 



1/ T, °K' 1 


Figure 1.- Carbon filament-oxygen reaction at p_ = 10 ^ torr. 

°2 

(After Blyholder, Binford, and Eyring (ref. 2).) 
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(Determining reaction rates at low pressures has the virtue that the low pressure effec- 
tively eliminates transport limitations by increasing the mean free path of the gaseous 
reactant. For instance, at an oxygen pressure of 10“^ torr and a temperature of 2000° K, 
the mean free path of molecular oxygen is on the order of 8 cm.) A most interesting 
aspect of the data in figure 1 is that, without exception, maxima exist in the various rate 
curves. In addition, the data of Meyer (ref. 3) and Strickland-Constable (ref. 6) show that 
as the temperature is increased beyond the rate maximum, a minimum is reached after 
which the rate increases again. It is unlikely that sublimation is responsible for this 
second increase in rate, since the theoretical free vaporization rate (i.e., the theoretical 
maximum vaporization rate) of carbon is only 4.6 X 10"^ g/cm 2 sec at 2000° K. (See 
ref. 34.) Quantitative differences in the data presented in figure 1 can be ascribed to 
differences in structural and chemical (i.e., impurities) properties of the various carbons 
employed (refs. 2 and 35) in addition to possible differences in experimental technique. 

The plot of the kinetic expression developed for these data by Blyholder et al. is shown in 
figure 1 as the dashed line. It may be given as 

rfa = 2.58 X lO^kjfpQ + 3.4 X 10" 10 k 2 (l - f) (3) 

2 

where 

k x = 1.6 x io 8 e" 3000 °/ RT 

k 2 = 2.2 x io 9 e" 8000 °/ RT 

k 3 = 4.2 x 10 24 e _17500 °/ RT 

f _ 2k 2 
2k 2 + k 3 

The theoretical basis on which this equation was developed is that there are two types of 
reaction sites which are thermally convertible one into the other according to a first- 
order rate expression. 

Another high-temperature kinetic expression is that of Strickland-Constable and 
coworkers (refs. 8 to 10) and is based on their data on the reaction rate of pyrolytic 
graphite. Experimentally, pyrolytic-graphite rods 4 mm on a face by 2 mm thick were 
electrically heated, and an unheated jet of oxygen or of oxygen-nitrogen mixtures was 
impinged on their surfaces at high velocities in the (impact) pressure range 0.20 to 
0.23 atm oxygen. Rate maxima were observed at about 2173° K. Similar experiments 
conducted on lampblack carbon and reactor -grade graphite also showed rate maxima but 
at absolute rates about nine times faster than those for the pyrolytic graphite. In view of 


6 



this result, the very compelling argument was advanced (ref. 8) that since oxygen can 
obviously be transported to the pyrolytic -graphite surface nine times faster than the rate 
required to maintain its observed reaction rate, its reaction rate clearly cannot be trans- 
port limited. This conclusion was further substantiated by an independent experiment 
(ref. 9) in which the reaction rate of the pyrolytic graphite was studied as a function of 
jet velocity at the temperature of the rate maximum. At low jet velocities, the reaction 
rate increased with velocity; but above 2500 cm/sec, the rate leveled off, no longer 
increasing with velocity. However, at the same experimental conditions, the rate of reac- 
tion of a reactor -grade graphite continued to increase with increasing velocity. Hence, 
it could be safely concluded that the observed reaction rate of the pyrolytic graphite 
represented the true kinetic rate and was not influenced by transport limitations. A 
kinetic expression was developed in reference 8 for the pyrolytic -graphite data on the 
assumption that CO is the primary product of the reaction and that there are two types 
of active sites on the surface which undergo thermal inter convertibility. This kinetic 
expression may be given as* 

^ _ 12k A k BPo 2 2 12k T k BPo 2 

(* + k zP0 2 )( k T + k BP0 2 ) k T + k BPo 2 

where 


k A = 20e" 3000 °/ RT 

k B = 4.46 x 10“ 3 e" 15 290//RT 
to* 

k T =1.51x«^e- 9 ' 700 °/ RT 


k z = 21.3e 410 °/ RT 


This expression was further shown to be in good qualitative agreement with the earlier 
carbon-filament data of Strickland-Constable (ref. 6) at oxygen pressures of 2.5 x 10“5 atm 
and of Duval (ref. 4) at oxygen pressures of 10 atm. 

Ong, in an impressive theoretical analysis (ref. 11), developed yet another kinetic 
expression valid at high temperatures. He applied the transition -state theory of kinetics 
to a deduced set of controlling kinetic steps, calculating the required transition energies 
by statistical thermodynamics, and evaluating several undetermined constants in his final 
expression from the data of Okada and Ikegawa (ref. 12). Ong’s theoretical development 

*The expression originally given by Nagle and Strickland-Constable (ref. 8) con- 
tained a typographical error which was subsequently corrected in references 9 and 10. 
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(ref, 11) differs considerably from those of Blyholder et al. (ref. 2) and Strickland- 
Constable and coworkers (refs. 8 to 10) in that he assumed only one type of active site 
on the carbon surface. Another important feature of Ong's development is that he con- 
sidered two product species from the carbon-oxygen reaction, CO 2 at low temperatures 
and CO at higher temperatures. Also, a rate maximum is predicted by Ong's kinetic 
expression, which he attributes to a change in the chemical potential of vibration of 
molecular oxygen. Ong's rate expression may be given as 


d. = F^n.e- 2280(, / RT 


; i p o 2 


1 + ( k iVo 2 ) + k i p o ; 


2F™e-51400/RT 


( k l k 3 p 0 2 ) : 

( k l k 3»02) 1A 


+ k iPr 


where 4 


F C02 = 2 - 95 

2F CO = 9.5 x 10 3 

kl = e -(- 78410 - p l)/ RT 

k 3 = e-(- 5760 - p 3)/ R T 

Mj = RT^to ^ kT 


M3-RT& + ln(l -e“ 8 v/ T ) 


2 Table 5 of Ong's paper (ref. 11) lists F^q^, as 2.95 x 10 and F^q (instead of 
2F CO ) as 9.5 x 10 3 . These values are apparently in error when compared with figure 1 of 
his paper. The corrected values are shown above. In the expression for kg, Ong gives 
the constant in the exponent as -7990 cal/mole, however, the value of -5760 cal/mole has 
been chosen here because it gives better agreement with figures 4 and 5 of Ong's paper 
and also with the rate data of Okada and Ikegawa as shown in figure 1 of Ong's paper. 
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where 


m = Molecular mass of oxygen 
k = Boltzmann constant 
h = Planck constant 
a = 2, a symmetry number 

0 r = 2.07° K, characteristic temperature for rotation 
c o e = 3, degeneracy of the electronic ground state 
Q v = 2230° K, characteristic temperature for vibration 

Since several constants in equation (5) were determined from the data of Okada and 
Ikegawa (ref. 12), it may be of interest to discuss briefly their experiment. Cylindrical 
specimens of artificial graphite, 10 mm in diameter, were heated electrically in room- 
temperature air at 1 atm with oxygen supplied only by natural convection. Surface tem- 
peratures were between 973° K and 2273° K, and Okada and Ikegawa assumed a transport- 
limited reaction. Although no experimental proof is given to support this assumption, it 
seems reasonable because of the experimental conditions. In view of this uncertainty, it 
is peculiar that Ong should have chosen these data for determining his kinetic constants. 
Nevertheless, it is interesting to note that whether or not the data are truly transport 
influenced, a very good representation of them can be obtained on the basis of theoretical 
considerations of reaction kinetics and statistical thermodynamics alone. It becomes 
apparent, then, that a nonlinear rate curve does not inevitably imply transport influence. 
This point has also been made by Rosner and coworker in references 18 and 19. 

In addition to the three kinetic expressions already presented, a fourth expression 
can be formulated by employing a suggestion of Essenhigh (ref. 13) in discussing rate data 
obtained by Golovina and Khaustovich (ref. 14). These data are for 15-mm-diameter 
spheres reacted in oxygen-nitrogen mixtures of air composition (presumably at 1-atm 
total pressure) at flow velocities of 60 cm/sec. Essenhigh deduced that at these high 
velocities, oxygen transport should not significantly influence the reaction rate at the 
experimental conditions of Golovina and Khaustovich. (This is partially confirmed by 
other data presented by Golovina and Khaustovich showing that the rate of the carbon- 
carbon dioxide reaction is independent of gas velocity in the 20 cm/sec to 80 cm/sec 
velocity range.) Essenhigh proposed that the rate data might be adequately represented 
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in terms of a simple adsorption -desorption mechanism with desorption controlling at low 
temperatures and adsorption controlling at high temperatures. The activation energy of 
desorption was assumed to be on the order of 40 000 cal /mole, and the activation energy 
for adsorption was reasoned to lie in the range 1000 to 5000 cal/mole. Although Essenhigh 
did not compare this proposed mechanism with Golovina and Khaustovich's data, this 
comparison can be made easily. By taking the activation energy for adsorption to be 
2000 cal/mole and by evaluating two undetermined constants from the data of Golovina 
and Khaustovich, a kinetic expression was developed 


. k l k 2Po 2 
m = k lP0 2 +k 2 

where 


( 6 ) 


kj = 5.48 x io- 3 e _200 °/ RT 
k 2 = 5.21 x 10 3 e' 4000 °/ RT 


Comparison of this expression with the data of Golovina and Khaustovich is shown in fig- 
ure 2. Although this expression fairly satisfactorily represents the data, it suffers from 
the oversimplification that it does not adequately account for the change in product from 
C0 2 to CO as temperature increases, nor does it predict a rate maximum at high temper- 
atures. What is important, however, is that a distinctly nonlinear rate curve (when plot- 
ted in an Arrhenius plane) can be adequately described by elementary considerations of 
heterogeneous kinetics without recourse to the concept of transport limitations. 



Figure 2.- Carbon sphere-air reaction. Data of Golovina 
and Khaustovich (ref. it) . 
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Equations (3) to (6) represent the only kinetic expressions which could be found or 
easily formulated from data in the literature. However, the literature also contains 
extensive additional data on the carbon-oxygen reaction at high temperatures further 
indicating distinctly nonlinear rate curves; these data are briefly discussed in appendix A. 


For purposes of comparing the predictions of the several kinetic expressions 
(eqs. (3) to (6)), mass-loss rates have been calculated at air pressures of 0.01 atm, 1 atm, 
and 100 atm and are plotted in figures 3 to 5. The temperature and pressure ranges of 
the supporting data for these kinetic expressions are listed in table 2. For comparison, 
the predictions of Scala (ref. 23) "fast" and "slow" kinetic expressions (eqs. (1) and (2)) 
and predictions of theoretical maximum stagnation -point transport -limited rates, based 
on the expression 


m D 


= 1.711 x 10" 2 



(7) 


from reference 23, for nose radii of 0.1, 1.0, 10, and 100 cm are also included in fig- 
ures 3 to 5. From these figures, it is seen that at temperatures above about 1700° K, 
all four high -temperature -rate expressions predict rates one order of magnitude or more 
below a linear extrapolation of the "slow" rate expression, and five orders of magnitude 
or more below a linear extrapolation of the "fast" rate expression, (it is of interest to 


TABLE 2.- KINETIC EXPRESSIONS 


Author 

Reference 

Equation 

Range of supporting kinetic data 

T, °K 

1 

PC>2’ atm 

Scala ("slow") 

23 

1 

698 to 1373 

See table 1 

Scala ("fast") 

23 

2 

698 to 1373 

See table 1 

Blyholder et al. 

2 

3 

1050 to 2300 

a 0. 0000132 

Strickland - 

8 to 10 

4 

1273 to 2673 

b 0.20 to 0.23 

Constable 





Ong 

11 

5 

973 to 2273 

0.21 

Essenhigh 

13 

6 

1073 to 2773 

c 0.21 

(extended) 






a Common pressure to which Blyholder et al. extrapolated data from various 
sources and at which they formulated their rate expression. 

■u 

D Impact pressure of the gas impinging on the carbon surface. The kinetic expres- 
sion was, however, also qualitatively substantiated at 2.5 x 10" 5 atm (refs. 8 and 10), and 
at 1 X 10"® atm (ref. 10). 

c The test-gas composition was equivalent to air and, although not stated, the pres- 
sure may be presumed to be atmospheric. 
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3000 TEMPERATURE, °K 

4000 2500 2000 1500 1250 1000 850 700 



17 T, 0 K _1 

Figure 3-- Comparison of carbon-oxidation-rate expressions 
at p a = 0.01 atm. 

note that the reaction rates predicted by the kinetic expressions in table 1, when com- 
pared with Scala's "fast" and "slow" rate bounds, are much closer to the "slow" bound 
with the one exception of the high rates predicted by the expression of Vulis. (For an 
explanation of why the rates predicted by Vulis are so high, see footnote (e) of table l] 

It is reassuring, then, that the present high -temperature -rate expressions form a more 
natural continuation of the "slow" rate bound than they do of the "fast" rate bound.) 

At pressures of 0.01 atm, none of the four high-temperature-rate expressions pre- 
dict a reaction rate fast enough to reach the maximum transport limits, even for a nose 
radius as large as 100 cm. At the higher pressures of 1 atm and 100 atm, and with the 

O 

exception of the curve of Blyholder et al., these rate expressions predict reaction rates 
below 2500° K which equal or exceed the transport limits only for relatively large nose 
radii. This result clearly indicates that for a wide variety of combinations of pressure, 
temperature, and nose radii, the maximum attainable rate will be determined by chemical 
kinetics and not by transport limitations. This statement is not to be interpreted to imply 

3 The reason the curve of Blyholder et al. (ref. 2) differs so noticeably from the 
others may be because Blyholder et al. developed their kinetic expression from data 
obtained at pressures far below the range of pressures in figures 3 to 5. 
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Figure 4 .- Comparison of carbon-oxidation-rate expressions 

at p a = 1 atm. 


that the carbon-oxygen reaction at high temperatures is never influenced or limited by 
oxygen transport. Certainly, at sufficiently low flow rates the reaction may well become 
limited by transport. But what is indicated is that at high temperatures the kinetics of 
the reaction are much slower than has generally been suspected on the basis of linear 
extrapolation of Scala's (ref. 23) "slow" and "fast" rate expressions; therefore, it cannot 
be assumed that a transport limitation inevitably must be reached. Further, to determine 
experimentally whether a reaction rate is chemically or transport controlled, the effect 
of environmental conditions of temperature, pressure, and flow velocity in relation to 
specimen size (and also carbon type) on the reaction rate must be considered; certainly 
no definitive statement predicting transport control can be made merely on the basis of 
the shape of the rate curve as a function of temperature. 

The kinetic expression which appears most reliable for estimating the oxidation rate 
of pyrolytic graphite at high temperatures is that of Strickland-Constable and coworkers 
(eq. (4)). For estimating the oxidation rates of other less graphitic carbons, the case is 
not so clear. The kinetic expression of Blyholder et al. (eq. (3)) is not as appropriate 
for use at normal pressures as the remaining expressions (eqs. (5) and (6)) because it is 
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Figure 5-- Comparison of carbon-oxidation-rate expressions 
at p a = 100 atm. 

based on data all obtained at oxygen pressures below 2 x 10“^ torr (some pressures were 
as low as 5 x 10"5 torr). These data were extrapolated in reference 2 to the common 
pressure of 10“^ torr in order to formulate a kinetic expression. No test of the resulting 
kinetic expression was made at other pressures. The kinetic expression developed by 
Ong for artificial graphite (eq. (5)), although particularly appealing because of its strong 
theoretical foundation, suffers from possible influences of transport interference because 
of its dependence on the rate data of Okada and Ikegawa for the evaluation of several 
constants. The kinetic expression given by equation (6), which is based on the data of 
Golovina and Khaustovich for an electrode carbon, although free of transport interference, 
suffers from certain theoretical oversimplifications as mentioned earlier. Nevertheless, 
its use appears preferred because it both fairly adequately describes the data and is 
simple to use. 


CONCLUDING REMARKS 

A large body of experimental data exists, convincingly demonstrating that the oxi- 
dation rate of carbon passes through a maximum with increasing temperature. Such a 
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maximum has been observed for such diverse carbons as pyrolytic graphite, artificial 
graphite, and vitreous carbon. At temperatures above the maximum, a rate decrease 
followed by a subsequent rate increase has often been observed. Thus, it is clear that 
the mere existence of a nonlinear rate curve on an Arrhenius plot is insufficient proof of 
transport influence in a reaction; in order to determine if a reaction is truly kinetically 
or transport controlled, effects of temperature, pressure, specimen size, flow velocity, 
and carbon type on the reaction rate must be considered. 

Four kinetic expressions based on kinetic data obtained ?t high temperatures have 
been presented and compared at air pressures of 0.01 atm, 1 atm, and 100 atm. Above 
about 1700° K, all these expressions predict mass -loss rates one or more orders of 
magnitude below that predicted by a linear extrapolation of Scala's "slow" kinetic expres- 
sion (Tech. Inform. Ser. No. R62SD72). At the pressures of these comparisons, the 
kinetic expression of Blyholder et al. (J. Phys. Chem., vol. 62, no. 3, 1958) is felt to be 
of doubtful reliability because it was formulated from data obtained at considerably lower 
pressures. Comparison of the other three rate predictions with theoretical maximum 
stagnation -point transport-limited rates of Scala (Tech. Inform. Ser. No. R62SD72) shows 
that only for relatively large nose radii will transport limitations be encountered below 
2500° K; that is, in many practical situations, the maximum attainable rate will be pre- 
dicted by kinetic considerations. Obviously, however, at sufficiently low flow rates the 
reaction may well be transport limited. 

For estimating the oxidation rate of pyrolytic graphite at high temperatures and 
normal pressures, the rate expression of Strickland-Constable and coworkers (Carbon, 
vol. 1, no. 3, 1964; and Second Conference on Industrial Carbon and Graphite, Soc, Chem. 
Ind. (London), 1966) appears tc be preferred: 

12k A k B p O ? 2 12k T k B r 0 

m = r-7 r- + : — 

l 1 + k z p 0 2 )( k T + k B p 0 2 ) T + k B p 0 2 


where 


kA = 20e -29 ppp / R,p 

k B = 4.46 x io- 3 e - 1520 °/ RT 
iCT* 

k T =1.51x|^e- 97000 / RT 
k z = 21.3e 410 °/ RT 
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m = Mass-loss or reaction rate, g/cra^sec 


= Oxygen pressure, atm 

For other less graphitic carbons, the rate expression developed in the present 
paper based on the data of Golovina and Khaustovich (Eighth Symposium (International) 
on Combustion, Williams & Wilkins Co., 1962) appears preferred: 

k l k 2 p 0 2 
m = - 

k l p 0 2 +k 2 


where 


kj = 5.48 x io- 3 e' 200 °/ RT 
k 2 = 5.21 x io 3 e _4000 °/ RT 
rh = Mass-loss or reaction rate, g/cm 2 sec 



= Oxygen pressure, atm 


As a final point, it is recognized that certain contradictory aspects exist in the data 
discussed, for example, the location of the rate maximum, the magnitude of the reaction 
rate, and effects of thermal accommodation, and furthermore, that a general uncertainty 
exists regarding the true kinetic mechanism. Certainly, an understanding of the carbon- 
oxygen reaction at high (or even low) temperatures is far from complete. Particular 
areas requiring study include the effects of inert diluents, the determination of pressure 
dependence over a wide range of pressure (particularly at higher pressures), the role of 
thermal accommodation and reaction hysteresis, as well as the problem of why the reac- 
tion rates of diverse carbons often differ so markedly. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Hampton, Va., May 5, 1971. 
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APPENDIX A 


ADDITIONAL REACTION-RATE DATA AT HIGH TEMPERATURES 

In addition to the kinetic expressions and data already discussed, results of a num- 
ber of other studies of the carbon-oxygen reaction at high temperatures have also been 
published. (See refs. 15 to 22.) These studies lend further support to the argument that 
reaction rate, and not oxygen transport, is often the rate -controlling mechanism. Fur- 
thermore, a rate maximum was observed in all of these studies. One of the more impres- 
sive of these studies is that of Lewis, Floyd, and Cowlard (ref. 15) who approached their 
study with an expressed concern that perhaps observed rate maxima were not truly asso- 
ciated with the reaction itself but were produced by experimental techniques. Of partic- 
ular concern to them in this regard was the possible effect of thermal accommodation 
(i.e., cold gas and hot surface). Results of their experimental study did show that ther- 
mal accommodation had an effect on reaction rate (higher gas temperatures resulted in 
higher reaction rates at the same specimen surface temperature), but the shape of the 
reaction-rate curve did not change significantly; that is, rate maxima were still observed. 
Although the study of thermal accommodation effects was conducted on the carbon-carbon 
dioxide reaction, it may be expected that the carbon-oxygen reaction would behave quali- 
tatively similarly. In the oxidation of vitreous carbon in a 0.1 O 2 and 0.9 argon mixture 
at 2 atm, a rate maximum occurred at 2050° K and a minimum at 2750° K. Other studies 
of the oxidation rates of pyrolytic graphite, vitreous carbon, and artificial graphite in CO 2 
demonstrated that both the location and the magnitude of the maximum are dependent on 
carbon type. 

Other interesting reaction-rate phenomena have been reported by Rosner and 
Allendorf (refs. 16 to 20) who studied the reaction of carbon with both molecular and 
atomic oxygen in the temperature range 1100° K to 2100° K at oxygen pressures of 
10"2 torr. In addition to the low pressure of their tests which minimizes transport 
influence, transport effects were ruled out on the basis that flow velocity had no effect 
on reaction rate. (See refs. 17 and 20.) Their results, expressed in terms of oxidation 
probability, ^ showed that (1) a rate maximum exists for both pyrolytic and artificial 
graphite in both molecular and atomic oxygen, (2) atomic oxygen is more effective in 
removing carbon atoms than is molecular oxygen, (3) artificial graphite is more reactive 
than the c-face of pyrolytic graphite, and (4) the highest oxidation probability in any of 


^Defined as the ratio of the number of carbon atoms removed from the surface to 
the total number of reactant collisions with the surface. In two early papers (refs. 16 
and 17), the magnitude of the oxidation -rate probability was erroneously reported too 
high. This was corrected in subsequent papers. (See refs. 19 and 20.) 
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APPENDIX A - Concluded 


the tests was about one -half. This last point is of particular importance because it says 
that, even at the high temperatures of the study, a maximum of only one-half as many 
carbon atoms are being removed from the surface as the number of oxygen atoms striking 
the surface. When the oxidation probabilities of Rosner and Allendorf (ref. 20) are 
converted to conventional reaction rates and compared with the reaction-rate data of 
Blyholder et al. (ref. 2) shown in figure 1, their pyrolytic -graphite data closely follow 
the data of Sihvonen, and their artificial -graphite data describe a curve similar in shape, 
but higher than that of Eucken. 

A particularly interesting hysteresis effect in the oxidation rate of carbon has been 
reported by Duval (ref. 21) who reacted graphitic and nongraphitic carbon filaments in 
oxygen in the pressure range 5 x 10“® to 3 x 10"^ torr. In addition to the usual rate max- 
ima, he observed that a carbon, having once reached a constant reaction rate at a given 
temperature, will, when reacted at a new temperature, react at an initial rate which is 
higher or lower than the final steady-state rate at the new temperature, depending on 
whether the new temperature is higher or lower than the previous reaction temperature; 
that is, the initial rate of reaction at a given temperature was found to depend on the con- 
ditions of reaction immediately preceding. This result may be of great importance when 
dealing with rate transients due to changes in surface temperature. Another report of a 
rate maximum in the oxidation rate of carbon was given by Audubert and Busso (ref. 22) 
in their study of the reaction rate of electrically heated carbon filaments with air at 
20 torr. The rate maximum occurred at about 1825° K. 
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